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Abstract 

When initially-independent subsystems are made to contact, coherence can develop due to inter- 
action between them. We exemplify and demonstrate this paradigm through several scenarios of 
two initially-independent Bose-Einstein condensates which are allowed to collide. The build-up of 
coherence depends strongly on time, interaction strength and other parameters of each condensate. 
Implications are discussed. 
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Coherence is a fundamental property in quantum mechanics, and is the basis for both 
fundamental and practical research fields in modern physics [l, 2]. Often one speaks of the 
loss of coherence when a system is coupled to an environment, or a bath. In this context 
and for physical applications relying on coherence, one is also interested in how to protect 
coherence of quantum systems [3(. 

The issue we address in the present work is a complementary paradigm, namely build-up of 
coherence. We show that when initially-independent subsystems are made to contact, mutual 
coherence can develop due to interaction between them. We exemplify and demonstrate this 
paradigm through several scenarios of two initially-independent Bose-Einstein condensates 
(BECs) which are allowed to collide. It has been suspected before that coherence builds- 
up in the collision between two initially- independent BECs [4j. Here, for the first time, 
we show explicitly by solving the time- dependent many-body Schrodinger equation and 
resorting to the reduced one-particle density matrix that coherence can build-up and how 
it happens. The build-up of coherence depends strongly on time, interaction strength and 
other parameters of each BEC. 

We consider two initially-independent BECs, an A BEC with Na atoms and a B BEC 
with Nb atoms. In the simplest case, each of the two BECs is a weakly interacting Bose 
gas associated with an orbital, 4>a{ x ) for the A BEC and 4>b{x) for the B BEC. Being 
independent means that the two BECs do not overlap in space, J dx\(j)A{x)\ 2 \(j)B{x)\ 2 = 0. 
Nowadays, independent BECs can readily be realized if the two BECs are held in two 
different traps, or in a double-well potential with a large barrier between the wells. When 
the A and B BECs are comprised of the same kind of bosons, N = Na + Nb, the quantum 
state of the whole system reads ^ = S4>a{x\) ■ ■ ■ 4>a(xn a )4 ) b(xn a +i) ■ ■ ■ 'Pb (x n a + Nb ) > where 
S is the symmetrization operator. The many-boson state \I/ thus describes a fragmented 
BEC (5, 6, 7, 8|. 

At time t = we remove the traps and allow the whole system to evolve in free space 
under the many-body Hamiltonian H = + J2j<k W(xj — Xk). Here, T(x) is 

the kinetic energy operator and W(x — x') the inter-particle interaction. The system is 
taken to be one-dimensional. The many-body state of the system evolves according to 
the time-dependent Schrodinger equation Hty(t) = i^p-- We solve the Schrodinger equa- 
tion with the recently developed multi-configurational time-dependent Hartree for bosons 
(MCTDHB) [9]. In the MCTDHB(M) the time-dependent many-boson wavefunction #(t) 
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is written as a linear combination of all possible permanents resulting by distributing iV 
bosons over M time-dependent orthogonal orbitals. The MCTDHB(M) wavefunction reads 

ty(Xi,X 2 , ...,X N ,t) = J^rn Cm(t)S(j)i(Xi,t) ■ ■ ■ <fa{x mi+1 ,t) ■ • • (j) 3 (x mi+m2+ i, t) ■ ■ ■ <f) M {xN, t) , 

where in = (m l5 m 2 , • • ■ m^/) collects the occupation numbers. Within this theory a 
quantitative description of the time evolution of bosonic systems is achieved by opti- 
mizing all M orbitals used to construct the many-body expansion and the expansion 
coefficients themselves at each point in time utilizing a standard time-dependent varia- 
tional principle 9|]. To analyze the system's evolution we will also resort to the respec- 
tive time-dependent multi-orbital mean-field (TDMF) dynamics 1(3] of the two initially- 



independent BECs ll|, |l2j. In TDMF(M) theory the many-boson wavefunction has the 



form ...,x N ,t) = S<pi(xi,t) ■ ■ ■ <p 2 (x mi+1 , t) ■ ■ ■ <p 3 (x mi+m2+1 , t) ■ ■ ■ <fi M (x N ,t) : i.e., the 

time-dependent many-body wavefunction is comprised of one permanent only and thus main- 
tains the initial occupation numbers m in time. Obviously, the many-body dynamics goes 
much beyond the TDMF one. 

For a system comprised of identical bosons the most basic quantity quantify- 
ing coherence in the system is the reduced one-body density matrix p(x,x',t) = 
N J dx 2 ■ ■ ■ j cIxn\^/*(x', x 2 , ■ ■ ■ , xn, t)^(x, x 2 , ■ ■ ■ , x^, t)\. The reduced one-body density 
matrix has the familiar spectral resolution p(x,x';t) = ^2jnj(t)ijj*(x' ,t)ipj(x,t) where nj(t) 
are the natural occupation numbers and ipj (x, t) are the natural orbitals. For the scenarios 
studied here there are two principal natural orbitals occupied. Furthermore, we will dis- 
cuss the density p(x,t) = p(x,x = x',t), namely the diagonal part of the reduced one-body 
density matrix. 

We consider two-initially independent BECs with Na = Nb = 500 bosons each. This 
system is prepared as the ground-state of a symmetric double-well potential formed by 
bisecting an harmonic potential with a Gaussian-shaped barrier. Here and hereafter we 
work in dimensionless units which are arrived at by choosing a convenient length scale L 
(say the initial distance between the BECs) and dividing the Hamiltonian H by ^2 where m 
is the boson mass. The kinetic energy operator then reads T = —2^25 an d the double- well 
is V(x) = 0.05x 2 + 50/ V2n exp[—x 2 /8]. As the particle-particle interaction we employ the 
standard delta-function potential, W(x — x') = \q5(x—x'), where the transverse confinement 
is accounted for in Ao [13|. Here Ao = 0.1 

In Fig. [T] we present snapshots of the many-body density p(x, t). Also shown is the mean- 
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field density p MF (x,t) computed by the TDMF approach. Both densities show interferences 
of the two initially-independent BECs, a subject which has drawn much attention recently 
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. At short times, up to about t = 3.0, the two densities and thus the interference 



patterns coincide. Thereafter, we see that differences start to appear between the two 
densities. It has been shown recently that interaction leads to interferences in the density 



of two initially-independent BECs 
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121 ] . The interaction between the two expanding 



12|. This 



BECs leads to "interaction-assisted self-interference" mechanism between them 
mechanism is seen in the mean- field density p MF (x, t) of Fig. [TJ The differences between the 
many-body and mean-field densities seen in Fig. [T]thus indicate that more than "interaction- 
assisted self-interference" is happening in the system. 

To understand and quantify the differences between the many-body and mean-field den- 
sities we begin by plotting in Fig. [2] the natural occupation numbers Uj(t) of the many-body 
solution. We recall that in the mean-field dynamics the natural occupation numbers do not 
change in time. The natural occupations numbers nj(t) stay "flat" at their initial values 
up to about t = 4. Note that the two initially-independent BECs have already met at 
that time, see Fig. [U From about t = 4 on, the natural occupation numbers change from 
their initial values and coherence develops in the system. We see that coherence develops 
in an oscillatory manner and not monotonously. At about t = 7 the coherence assumes its 
maximal value for the system under investigation. It is instructive to employ the quantity 
C(t) = 100 x "^""(o^ which indicates the fraction in percents at time t of the maximally 
possible coherence in the system. By definition, C = at t = 0. In our system of two 
initially- dependent BECs with the same number of bosons C(t) takes a simplified and ap- 
pealing form, C{t) = 100 x , measuring the total coherence in the system. For two 
initially-independent BECs C — 0, and if the whole system becomes condensed C = 100%. 
From Fig. [2] we see that in our example the maximal coherence developed (t ~ 7) in the 
whole system is about C = 34%. 

The result depicted in Fig. [2] clearly shows that coherence builds-up in the system as 
a consequence of the interaction between the initially-independent BECs. The coherence 
changes with time in a nontrivial oscillatory manner. We have also examined the build-up 
of coherence in a smaller system of N = 100 bosons. We considered two initially-independent 
BECs with Na = N B = 50 bosons prepared in the ground-state of the above double-well 
potential V(x). The interaction strength Aq is tuned such that the factor Xq(N — 1) is the 
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same as for the above-studied case with N = 1000 atoms. The build-up of coherence is 
clearly observed in this smaller system as seen in Fig. [2j It is interesting to note that the 
build-up of coherence of the two systems follows a similar time-dependent pattern. 

It is instructive to return to the density in Fig. [T] and analyze how the build-up of 
coherence influences the density p(x,t). As there are two principal natural occupation 
numbers the density is expressed as 



The initial conditions of the two initially-independent BECs are: ni(0) = Na, ipi(x,0) = 
4>a{x) and 712(0) = N B , ?p2{x, 0) = <$b{x). With time, the interaction between the bosons 
modifies the natural orbitals ijji(x,t),ijj2{x,t) as well as the natural occupation numbers 
n 1 (t),n 2 (t). Similarly, within the TDMF approach [lOj the density is composed of two 
natural orbitals as follows 



The initial conditions are the same, (/>a(x, t) = 4>a{ x ) an d ^(x, 0) = but the occu- 

pation numbers remain fixed-in-time. 

To show how many-body physics and the build-up of coherence manifest themselves 
in the time-dependent density, it is deductive to relate the two densities p(x, t) and 
p MF (x,t). For this, we expand the many-body natural orbitals i/ji(x,t),i/j2{x,t) which are, 
of course, normalized and orthogonal to one another with the help of the TDMF orbitals 
cj) A (x,t),(j) B (x,t). Specifically, we write ipi(x,t) = ^1 - \A 1B (t) 2 \ - \A 1Q (t)\ 2 (f)A(x, t) + 



A 1B (t)M^t) + A 1Q (t)0 1Q (x,t) and fo{x,t) = ^1 - \A 2A {t)\ 2 - |A 2Q (t)| 2 B (x, t) + 



A2A(t)(pA(x, t) + A 2 Q(t)(j)2Q(x, t). The first two terms in each expansion are the TDMF 
orbitals which are themselves orthogonal to one another. The remaining term of each ex- 
pansion 0iq(x, t), 4>2q(x, t) belongs to the subspace orthogonal to TDMF orbitals. The 
functions 0iq(x, t) and <p2Q{x,t) are in general non-orthogonal. Similarly, we write for the 
time-dependent natural occupations rii(t) = Na + An(t) and n 2 (t) = N B — An{t). Thus, by 
construction, at t = all quantities denoted with A are equal to zero; Their non-zero values 
obtained with time signify many-body facets of the interaction in the system. 

It is straightforward to insert the above expansions for the natural orbitals and occupa- 
tion numbers into Eq. ([1]). The expression obtained is lengthy and at first sight not very 



p(x, t) = ni (t) |^i (x, t)\ 2 + n 2 (t) t)\ 2 . 



(1) 



p MF (x,t) =N A \M^t)\ 2 + N B \Mx,t)\ 2 - 



(2) 
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informative. More informative is to analyze p(x, t) to first order in the quantities denoted 
in the expansions by A. The result takes on the following form: 

p(x,t) = p MF (x,t) + An(t) [|<MM)| 2 -|<MM)| 2 ] 

+ 2Re{N A [A 1B <l> m A (x,t)<f>B(x,t)]+NB [A 2A <f) B (x,t)<f) A (x,t)}} 

+ 2Re {N A [A 1Q( f> A {x, t)<f> 1Q {x, t)\ + N B [A 2Q <t>* B {x, t)(j> 2Q {x, t)}} . (3) 

With its four terms, Eq. ([3]) offers an appealing and explicit mechanism for the many-body 
effects of interaction on the density of two initially-independent BECs. The first term is 
just the mean-field density p MF (x,t). As discussed above and shown in Fig. [TJ at short 
times this is the only (visible) contribution to the density. As time progresses corrections 
to the mean- field quantities start to mount, also see Fig. [TJ resulting in three additional 
contributions to p(x,t). The second term in ([3]) is due to the build-up of coherence in the 
system, and describes the flow of particles from one BEC to another without changing the 
TDMF orbitals themselves. The third term in ([3]) is an interference term and describes 
interference between the two BECs <f) A (x,t) and (pB{x,t). It originates from the change 
of the many-body natural orbitals within the TDMF subspace {<p A (x, t), t)}, without 
changing the TDMF occupation numbers. The last, fourth term is another interference term. 
The many-body natural orbitals start to change form beyond the linear-combinations of 
TDMF orbitals. The fourth term thus describes the interference of the A BEC [4> A (x, t)} with 
4>iQ(x,t) and, separately, the interference of the B BEC [(j) B (x,t)] with cj) 2 Q(x,t). Adding up 
the above four contributions, the time- dependent density of two initially-independent BECs 
reflects combined interference and build-up-of-coherence mechanisms. Finally, we remark 
that going beyond first order in the A terms, see Eq. ([3]), there are more interference terms 
and terms resulting from the build-up of coherence. Of course, the many-body dynamics 
computed with the MCTDHB approach and depicted in the figures are to all orders in A. 
The above analysis in A is for the sake of physical interpretation only. 

Next, we would like to study more on the development of coherence in a system with two 
initially-independent BECs. We have seen that coherence develops in time and the question 
we would like to address is whether we can influence this process. We make the following 
'experiment'. We release the two BECs and give the A BEC velocity v and the B BEC 
the opposite velocity —v. Positive v means that the BECs move towards each other and 
negative v the reverse. The relative velocity of the two BECs is of course 2 x v. Examples are 
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collected in Fig.[3]where the natural occupation numbers rij(t) of the many-body solution are 
plotted. The faster the two BECs move towards each other, the quicker coherence initially 
develops (see the maximal rii(t) at about t = 7, 5,3 for velocities v = 0, 1,2, respectively). 
This behavior can be anticipated by the kinematics: the faster the BECs move towards each 
other, the sooner they come in contact which allows coherence to build-up. Side-by-side, we 
observe that the coherence oscillates and with increasing positive velocity v, the height of 
the first oscillation decreases. An interesting point seen in Fig. [3] for positive v is that this 
decrease is approximately linear in v. Let us make another 'experiment', in which the two 
initially-independent BECs move away from one another. The natural occupation numbers 
are plotted for v = —2 in Fig. We first see that coherence develops slower than for positive 
v, in conjunction with the above kinematic analysis. However, its maximum (C ~ 42%) is 
even higher than the maximal amount of coherence developed with initially-independent 
BECs at rest (v = 0). To explain this, at first sight counterintuitive result, we recall that 
the build-up of coherence between initially-independent BECs implies that more bosons 
share the same natural orbital ipi(x, t). How to 'help' initially-independent bosons share the 
same natural orbital? Bosons coming from initially-independent BECs at rest have positive 
relative velocity due the free expansion of the two BECs. These bosons are more 'distinct' 
from one another and therefore less coherence can build-up. Sending the two BECs with 
negative relative velocity compensates partly for this expansion of the two BECs. The now 
slower bosons are more 'similar' to one another and therefore higher degree of coherence 
builds- up. 

Let us summarize. Coherence can develop between two initially-independent BECs. How 
it develops in time when the two BECs are released from their trap strongly depends on 
different parameters, such as the inter-particle interaction and relative velocity between the 
BECs. We have considered the build-up of coherence between two initially-independent 
BECs comprised of indistinguishable particles. This motivates the consideration of other 
systems made of initially-independent subsystems where we also may anticipate build-up of 
coherence in time. For instance, the possible build-up of coherence between two subsystems, 
one coherent the other not, is an interesting case. In the context of quantum gases, the 



two subsystems can be a BEC and a Mott-insulator in an optical lattice [15|, [16| . Another 
appealing direction is when the initially-independent subsystems are made of different kinds 
of particles. Here, due to the distinguishablity of the particles, it would be instructive to 
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enquire on the possible build-up of coherence within each subsystem as well as in the system 
as a whole. The latter is measured by higher-order density matrices. 
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FIG. 1: (Color online) Free expansion and interferences of two-initially independent BECs each 
made of 500 bosons for Ao = 0.1. Shown as a function of time are snapshots of the density. On the 
multi-orbital mean-field level [p MF (x,t) - red curves] interferences result from "interaction-assisted 
self-interference" mechanism of the two BECs [ijj]. On the many-body level \p(x, t) - black curves] 
there is in addition build-up of coherence in the system. See text for more details. The quantities 
shown are dimensionless. 
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FIG. 2: (Color online) Build-up of coherence between two initially-independent BECs. Shown 
are the natural occupation numbers rij(t) as a function of time for the system in Fig. [TJ For 
comparison, the evolution of rij(t) for two initially-independent BECs with 50 bosons each and the 
same \q(N — 1) factor is shown. The quantities shown are dimensionless. 
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FIG. 3: Influence of relative velocity 2 x v on the build-up of coherence. Positive v indicates that 
the BECs move towards each other and negative v the reverse. Shown are the natural occupation 
numbers rij(t) as a function of time for two initially-independent BECs with 50 bosons each. The 
quantities shown are dimensionless. 
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